UNITU-THEP-2/1998 



FAU-TP3-98/1 



Solving the Gluon Dyson— Schwinger Equation 
in the Mandelstam Approximation 



A. Hauck^, L. von Smekal^, R. Alkofer^ 

'^Auf der Morgenstelle I4, Institut fur Theoretische Physik, Universitdt Tubingen, 

72076 Tubingen, Germany 

^Institut fiir Theoretische Physik III, Universitdt Erlangen-Niirnberg, Staudtstr. 7, 

91058 Erlangen, Germany 



Abstract 

Truncated Dyson-Schwinger equations represent finite subsets of the equations of 
motion for Green's functions. Solutions to these non-hnear integral equations can 
account for non-perturbative correlations. We describe the solution to the Dyson- 
Schwinger equation for the gluon propagator of Landau gauge QCD in the Mandel- 
stam approximation. This involves a combination of numerical and analytic meth- 
ods: An asymptotic infrared expansion of the solution is calculated recursively. In 
the ultraviolet, the problem reduces to an analytically solvable differential equation. 
The iterative solution is then obtained numerically by matching it to the analytic 
results at appropriate points. Matching point independence is obtained for suffi- 
ciently wide ranges. The solution is used to extract a non-perturbative /3-function. 
The scaling behavior is in good agreement with perturbative QCD. No further fixed 
point for positive values of the coupling is found which thus increases without bound 
in the infrared. The non-perurbative result implies an infrared singular quark in- 
teraction relating the scale A of the subtraction scheme to the string tension a. 

PACS Numbers: 02.30Rz 11.15.Tk 12.38.Aw 14.70.Dj 



Preprint submitted to Elsevier Preprint 



1 February 2008 



PROGRAM SUMMARY 



Title of program: mandelstam 
Catalogue identifier: 

Program obtainable from: CPC Program Library, Queen's University of Belfast, 
N. Ireland 

Computers: Workstation DEC Alpha 500 

Operating system under which the program has been tested: UNIX 

Programming language used: Fortran 90 

Memory required to execute with typical data: 200 kB 

No. of bits in a word: 32 

Peripherals used: standard output, disk 

No. of lines in distributed program, including test data, etc. : 347 
Keywords: 

Non-perturbative QCD, Dyson-Schwinger equations, gluon propagator. Lan- 
dau gauge, Mandelstam approximation, non-linear integral equations, infrared 
asymptotic series, constrained iterative solution. 

Nature of physical problem: 

An approach to describing non-pertubative correlations in field theories is to 
investigate their Dyson-Schwinger equations in suitable truncation schemes. 
Thereby one generally encounters non-linear integral equations which in many 
cases have infrared singular solutions imposing stability problems in the nu- 
merical procedures. 

Method of solution: 

The infrared singular part of the integral equation is studied analytically and 
transformed into constraints. Infrared subleading conributions are then ob- 
tained in form of an asymptotic series. This expansion of the solution can 
be calculated recursively and, together with the ultraviolet behavior obtained 
from a differential equation, it makes a numerical solution tractable. 

Restrictions on the complexity of the problem: 

So far all contributions from Faddeev-Popov ghosts are neglected. An ex- 
tended model based on analogous techniques is presented elsewhere [1,2]. 

Typical running time: Approximately 10 sec. 
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LONG WRITE-UP 



1 The physical problem 



1.1 Introduction 



The knowledge of the infrared behavior of the running coupling asifJ^"^) in 
QCD is crucial for an understanding of confinement. Should the theory have 
no further fixed point for > then the running coupling increases without 
bound in the infrared, and confinement might be realized by an absence of the 
cluster decomposition property in QCD. 

Due to the intrinsically non-perturbative nature of the problem, very little 
is known about the strong coupling as in the infrared. This is in contrast 
to its accurate knowledge in the ultraviolet, obtained via the perturbative 
calculation of the Callan-Symanzik f3 function. The logarithmic decrease of 
as signals asymptotic freedom and is verified experimentally in an impressive 
manner. 

Non-perturbative methods are required to study the strong coupling 05 in 
the infrared. One important framework for non-perturbative studies of QCD 
is provided by the infinite tower of its Dyson-Schwinger equations. These 
studies rely on specific truncation schemes of this tower. In the present paper 
we focus on an approximation scheme for the gluon Dyson-Schwinger equation 
in Landau gauge originally proposed by Mandelstam [1]. The corresponding 
equation will be referred to as Mandelstam's equation in the following. Two 
distinct approaches to its numerical solution are reported in the literature 
[2,3]. 

It was already pointed out by Mandelstam that for selfconsistent solutions to 
exist, certain conditions on the behavior of the gluon propagator in the infrared 
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have to be met [1]. An existence proof, a discussion of the singularity structure 
and an asymptotic expansion in the infrared for the solution to Mandelstam's 
original equation can be found in [2]. Taking care of various violations of gauge 
invariance in the approximation, one arrives at a slightly different equation [3]. 
We present an approach to the numerical solution of both equations, which 
combines methods of refs. [2] and [3]. The general arguments that follow apply 
to both versions of the equation. 

To remove the ultraviolet divergences in the equation for the gluon self-energy, 
wave-function renormalization is required, thus introducing an arbitrary scale 
ji. We demonstrate, how the rcnormalized equation can be cast in a renormal- 
ization group invariant form. We solve this equation numerically and deter- 
mine the running coupling from the solution via the renormalization condition. 
We show that the product of the coupling and gluon propagator, gD^^{k), 
does not acquire multiplicative renormalisation in the Mandelstam approxi- 
mation. This allows to identify a physical scale, the string tension, from the 
infrared singular result for this quantity. Consequently, our non-perturbative 
and renormalisation group invariant results yield a relation between the string 
tension and the only parameter in our calculation, the QCD scale A. This re- 
lation is in reasonable agreement with the respective phenomenological values. 
We obtain a non-perturbative f3 function and recover the scaling behavior of 
pcrturbativc QCD at the ultraviolet fixed point modulo small corrections due 
to the presence of ghosts in Landau gauge. For positive values of the coupling 
no further fixed point exists. The running coupling increases without bound 
in the infrared. 

The situation changes dramatically when ghosts arc included [4,5]. Then the 
running coupling approaches an infrared stable fixed point. As the correspond- 
ing calculations are very involved, however, we demonstrate the numerical 
method within the less complex case of solving Mandelstam's equation in this 
paper first. A generalization to the solution of the coupled gluon-ghost equa- 
tions will be presented elsewhere [6]. Nevertheless, the central idea of how to 
treat anticipated infrared singularities in non-linear integral equations is much 
more transparent and easier developed from the simplified equations discussed 
here. The present results allow furthermore for a detailed comparison of the 
two schemes. In particular, this allows to asses the infiuence and importance 
of ghosts in Landau gauge. 



1.2 Mandelstam's approximation 

Besides all elementary two-point functions, i.e., the quark, ghost and gluon 
propagators, the Dyson-Schwinger equation for the gluon propagator also in- 
volves the three- and four-point vertex functions which obey their own Dyson- 



4 



-1 



-1 



^^^^ g 



F(r(p 



Fig. 1. Diagrammatic representation of the gluon Dyson-Schwinger equation in 
Mandelstam's approximation. 

Schwinger equations. These equations involve successive higher n-point func- 
tions. Typical truncation schemes for this infinite tower rely on neglecting all 
higher but three-point functions expressing the latter in terms of propagators. 
This can be done with additional sources of information like the Slavnov- 
Taylor identities, which are entailed by gauge invariance. In the present study 
we neglect fermions considering a pure Yang-Mills theory, which is believed 
to reflect characteristic features of QCD. Studies by Brown and Pennington 
[7] indicate that the inclusion of quarks results in a suppression of the infrared 
part of the gluon propagator, however, no qualitative changes were observed. 
In the Mandelstam approximation ghosts are also neglected, because their con- 
tribution was anticipated to be small [1,3]. As a result, one obtains a simplified 
equation for the inverse gluon propagator in momentum spaceQ, 

D-^^^(fc) = D-,''^{k) + §^Jd\ rr (fc, q - k, -q) 

X D''\q)DP''{k - q)Tl^'"'{q, k - q, -k) , (1) 

where Dq and Fq are the bare gluon propagator and the bare three-gluon 
vertex, and F is the fully dressed vertex (we will use iVc = 3 colors). This 
together with a particularly simple form for the three-gluon vertex is the 
Mandelstam approximation [1]. 

The form for the three-gluon vertex is motivated by its Slavnov- Taylor iden- 
tity, which neglecting all ghost contributions in the covariant gauge reads: 

where p, q and k are incoming momenta and Z{k'^) is the invariant function 
of the dressed gluon propagator, 

Cnt^) = ^(*'"'-^)+?^ . (3) 



^ We work in Euclidean space with a positive-definite metric 5^"^ = diag(l, 1, 1, 1); 
color indices are suppressed. 
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The Slavnov-Taylor identity fixes the longitudinal part of the vertex. It can 
be derived along the general procedure of ref. [8]. In the simple version of Eq. 
(2) without ghost contributions the solution is: 

rr(P, k) = qy^'^ip - qy + A^{p\ q^^^ip + q)" 

+ q^S'^'^pq - p''q^){p - q)" + cycl. permut., (4) 



with 



(5) 



Assuming that Z{q^) is a slowly varying function one may write 

rr(P, k) = A^{p\ g^)rr(p, k) . (6) 

While this form for the full three-gluon vertex simplifies equation (1) even 
more than the use of another bare vertex, F = Fq, it is regarded superior to 
the latter since it accounts for some of the dressing of the vertex as it results 
from the Slavnov-Taylor identity. With F = F^ from (6) the Dyson-Schwinger 
equation can now be cast in the simple form 



D-''\k)^D-,''\k) + ^ ld\rr{k,q-k,-q) 

xD-^{q)D^{k-q)V^,''^{q,k-q-k) . (7) 

This equation is schematically depicted in Fig. 1. In solving it an additional 
condition is implemented in the Mandelstam approximation. Eq. (2) and its 
solution (4,5) for ^ entail that 

While imposing this additional condition seems consistent with the other as- 
sumptions in the Mandelstam approximation, it has to be emphasized that 
concluding (8) as a result of (2) rehes on neglecting all ghost contributions in 
covariant gauges. In fact, if ghosts are included in resolving the Slavnov-Taylor 
identity for the 3-gluon vertex, condition (8) has to be modified [4] . 

In a manifestly gauge invariant formulation the Dyson-Schwinger equation 
for the inverse gluon propagator in the covariant gauge would be transverse 
without further adjustments. The longitudinal part of the gluon propagator 
does not acquire dressing and cancels with the one of the bare propagator. 
This may be violated due to the neglect of ghosts, the violation of Slavnov- 
Taylor identities and also due to a regularization that does not preserve the 
residual local invariance, i.e., the invariance under transformations generated 
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by harmonic gauge functions = 0). The latter is the case for an 0(4) 

invariant Euchdean cutoff A, which we will use to regularize Eq. (7). In the 
present approximation, longitudinal terms can be eliminated by contracting 
Eq. (7) with the transversal projector, 

pi^-(k) = S'^'' - ^ . (9) 

After performing the angular integrations the equation for the gluon renor- 
malization function 2'(A;^) in Landau gauge = 0) then becomes [1] 



= 1 + 



Z(A;2) 167r2^ 




25 k ^ \ r-z / 9\ / 



fe2 



However, the above equation contains a quadratically ultraviolet divergent 
term. This term violates the masslessness condition (8) and has therefore been 
dismissed. As observed by Brown and Pennington [7], in general, quadratic 
ultraviolet divergences can occur only in the part of the inverse gluon prop- 
agator proportional to S'^". Therefore, that part cannot be unambiguously 
determined, it depends on the routing of the momenta. This is due to the var- 
ious violations of gauge invariancc mentioned above. The unambiguous term 
proportional to k^k''' can be obtained by contracting (7) with 



Rf^-{k) = 8^"" - 4^ . (11) 



In this case, upon angle integration instead of (10) one obtains [7] 



Z{k^) 16n^J k^ \2k^ 2 k^ 8 





+ 16^7^ 1 87 VJ ^ ^ 

It is clear that the assumptions made so far are drastic, and it is question- 
able how much of the original theory survived. Furthermore, inclusion of ghosts 
changes the conclusions on the infrared behavior of gluon propagator and the 
running coupling even qualitatively [4] . For the sake of a detailed understand- 
ing of what exactly can be attributed to ghosts, it is nevertheless necessary 
to thoroughly study the implications of the Mandelstam approximation, most 
remarkably on the running coupling, in order to compare to conclusions of the 
present and earlier work on the Dyson-Schwinger equations of Landau gauge 
QCD. Besides this, the major numerical concepts used in the present paper. 
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entirely devoted to the Mandelstam approximation, will apply to the coupled 
system of ghosts and gluons in Landau gauge with little changes. 

One might think that the particular problem problem with ghosts can be 
avoided using gauges such as the axial gauge, in which there are no ghosts in 
the first place. As for studies of the ghion Dyson-Schwinger equation in axial 
gauge [9-11] it is important to note that these rely on an ad hoc assumption 
on the tensor structure of the gluon propagator. The occurrence of a second 
independently possible term, which happens to vanish in perturbation theory, 
has been disregarded so far [12]. In fact, if the complete tensor structure of the 
gluon propagator in axial gauge is taken into account properly, one arrives at a 
coupled system of equations which is of similar complexity as the ghost-gluon 
system in the Landau gauge [13]. Progress is desirable in axial gauge as well, 
of course. A necessary prerequisite for this, however, is a proper treatment of 
the spurious infrared divergences which are well known to be present in axial 
gauge due to zero modes. This can be achieved by either introducing redundant 
degrees of freedom, i.e. ghosts, also in this gauge (by which it obviously looses 
its particular advantage) or by using a modified axial gauge [14], which is 
specially designed to account for those zero modes. Eventually, progress in 
more than one gauge will be the only reliable way to asses the influence of 
spurious gauge dependencies. 



1.3 Renormalization 



The logarithmic ultraviolet divergences in equations (10) and (12) can be re- 
moved by multiplicative renormalization. Introducing the renormalized gluon 

propagator Dji and the renormalized coupling g the renormalization constants 
Z3 and Zg are defined hj D = Z^Dji and Qq = ZgQ respectively. From the 
Dyson-Schwinger equation for the gluon propagator in Mandelstam approxi- 
mation (7), renormalized this way, we obtain 



(13) 



Here, /^^ := I[Zr] is a functional of Zjilk"^); for Eq. (12) Iz^ik"^) is given by 



1 




17 
2 A;2 



(14) 
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An analogous form follows for Mandelstam's original equation. After the sub- 
traction of the ultraviolet divergences, we use the limit A — ^ oo for the cutoff 
A implicitly. Confessing that we are dealing with renormalized quantities in 
the following, we dismiss the subscript R again. Below we will prove the iden- 
tity ZgZz — 1 for the Mandelstam approximation. Furthermore, we adopt a 
momentum subtraction scheme requiring the gluon self-energy to vanish at 
the renormalization point /x, i.e., 

z{^^') = l Z3 = i-^J^(/.^). (15) 

With this and asi/jt'^) = g'^ {fi^) / An Eq. (13) reads, 

(16) 

To determine the behavior of Z{k'^) for A;^ — > we first notice that 

(1) for a constant Z{k'^) the last term in (14) yields an infrared logarithm. 
Therefore, a solution to (16) cannot approach a constant as /c^ 0. Any 
constant term in the denominator on the r.h.s. of (16) has to vanish. This 
is a first constraint for a possible solution. 

(2) Similarly Z{k'^) ~ A;^ as A;^ — > is inconsistent since then the next to 
leading term yields a logarithmic contribution. 

(3) Now let us assume Z{k'^) ~ l/k"^. This form yields exclusively terms 
that violate the condition (8). Such terms have to be subtracted. Since 
the kernels of all integrals are linear in Z, this is achieved by simply 
subtracting a corresponding contribution from Z in the integrands. Thus 
defining F{k'^) by 

Z{k') = l,+ F{k') (17) 

with some constant 6, we further explore Ansatze for F, which is the 
remainder of Z in the integrals. 

To isolate the infrared singular term h/k"^ we rewrite equation (16) as: 

Z(.^)4 + ^(,^).^-^J-^, (18) 



z{e 



1 + 



An 
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where 



B 



8 47r 



7 J dq' 



dq 



2 F{<1') 



(19) 

(20) 

(21) 



As already stated in (1) above, a possible solution will have to obey A = 0. We 
will verify below that A = is not an independent constraint but an identity. 
In addition, we see that b = 1/B must be met, which will be imposed as a 
constraint on a possible solution. 

Introducing dimensionless variables and appropriately scaled gluon self-energy 
functions Z and F by 



47r 



X — 



b V Mi^') ■ 



and Z(x) 



X 



+ F{x) , 



(22) 



the equation under consideration can be rendered dimensionless and scale- 
independent. From equation (18) we obtain 



X 



1 + xF{x) 



A + Bx + 



dy (T ]p 17 y 9 X 7 x"^ 



n \ 



2x2 



(23) 



Here, 



B^bB 



Fjy) 

y2 



(24) 



from the constraint mentioned above. Since F{x) vanishes for a; — > as will be 
shown later, in this limit equation (23) shows explicitly that the x-independent 
constant 



T 47r , 
A = —A = . 

9 

With (24) equation (23) therefore yields 

x'^F{x) fdyfTy"^ 17 y 9 , _x 7 X 



(25) 



l + xF{x) 



X [2x^ 2 X S^^y Sy 



F{y) . (26) 
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Before we turn to the numerical solution of (26) with the constraint (24), some 
remarks need to be made. 



First note that the renormalization condition Z^jj?) — 1 implies that 

as{s) ^ At: Z\s) , (27) 

with s := \J? j)? and as(s) := q;5(//^). This is compatible with the definition 
of A which yields (L := AVA^) 



OLs{s) = 4:71 



2 

(28) 



i y 



The right hand side above is identical to the one of (23) with x ^ s, and thus 
(28) to (27). The fact that only /I = is consistent, i.e., the infrared behavior 
of the solution, entails that the running coupling is singular in the infrared. 

Furthermore, we observe that a solution to the system of equations (26) and 
(24) is a renormalization group invariant function. The argument for this is 
as follows: Assume that the product gDf^,^{k) does not acquire multiplicative 
renormalization. This means that Z oc gZ is a renormalization group invari- 
ant combination. For solutions of the form as given in (18) we conclude that 



as/{A'!r)b —: is a renormalization group invariant of mass dimension two. 
Therefore, the introduction of the dimensionless variable x = k'^/X^ in (22) 
does not introduce an implicit dependence on the renormalization scale /i in 
the explicitly scale-independent set of equations (26,24). This verifies the as- 
sumption, i.e., a solution Z{x) = gZ{k'^)/{4n) of (26,24) is a renormalization 
group invariant. For the renormalization constants this entails that ZgZ^ — 1 
in the Mandelstam approximation. This is to be compared to the Abelian ap- 
proximation to QCD, in which ZgZ^J^ = 1 or, equivalently, g'^D^y{k) is a renor- 
malization group invariant combination (one application of this is ref. [15]). 
Closer to the present case is the renormalization group improved one-loop re- 
sult in QCD. In perturbative QCD the power 6 of the coupling in the invariant 
combination g^D^y{k) depends on the number of quark-flavors Nf. Two ex- 
amples are 5 = 13/11 for Nf = and d = 1 for Nf = 3. Even though Nf = in 
Mandelstam approximation, it resembles the three flavor result at this point. 
However, in contrast to perturbation theory the identity ZgZ^ = 1{6 = 1) 
holds for all momentum scales // in Mandelstam approximation. 

We conclude the discussion of the renormalization with a comment on the 
choice of the non-perturbative subtraction scheme (15). Since gZ^k"^) is a 
dimensionless renormalization group invariant in Mandelstam approximation. 
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it is a function of the running coupling g{t, g), 

gZ{k^) = f{g{tk, g)), tk = ^ln k^fi^ . (29) 

Asymptotically, the momentum subtraction scheme is defined by f{x) — > x 
for X ^ 0. If the product gZ{k'^) is to have a physical meaning, e.g., as 
potential between static color sources, it should be independent under changes 
{g, fi) {g', fi') according to the renormalization group. Taken the idea of 
renormalization group invariance literally for arbitrary scales n', 

gZ{ii'^)^g'^g{\u{ii'/ii),g). (30) 

Then, f[x) = x, Wx, and, accordingly, Z{ii'^) = 1, W/i. This is the only 
physically sensible non-perturbative extension of the momentum subtraction 
scheme in the present context. Note that (29) with f{x) = x is identical to 

(27) with s ^ X, the renormalization scheme is thus equivalent to defining the 
non-perturbative running coupling by the product of the coupling constant 
and the gluon renormalization function in the Mandelstam approximation. 



2 Numerical methods 



2.1 Series expansion in the infrared and asymptotic behavior in the ultravio- 
lett 



A straightforward implementation of (26) fails due to delicate cancelations in 
the infrared. Therefore the expansion in the the infrared has to be pushed 
beyond leading order for Z{x). Making the Ansatz F{x) — ax'^ one sees that 
no solution exists for 7 < 1. For 7 > 1 eqn. (26) yields 

^2+7 = -Co^{-i)x^ - aCoo(7)a;'+'^ (31) 

with 

7 17 9 7 7 

^°'^^^^ = 2(7T^"2(^"8(7TI)+7"%^ ■ ^^^^ 

Obviously the coefficient 6*00(7) to vanish in order to permit a selfconsis- 
tent solution. This yields a cubic equation with a unique solution for 7 > 1: 

1 / 1099 \\ 13 



7o = -V229COS - arccos ^= a; 1.2705 . (33) 

' 9 \^3 V 229A/229yy 9 ^ ' 

An analogous condition was already found by Mandelstam for his original 
equation [1] . It is interesting to note that the value of 70 given here is numer- 
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ically very close to the one obtained by Mandelstam (70 = — 1 ~ 1.273), 
even though the coefficients in (12) and in (10) are completely different. 

Looking at the corresponding results given in [3] one observes that their in- 
frared leading terms for F(/c^) do not obey the constraint (33) for 70. We 
verified numerically that the expressions given in ref. [3] reproduce themselves 
under the integral equation in the momentum range given in the corresponding 
figure in this reference (Fig. 3 in [3]). However, trying to extend the numerical 
calculations further to the infrared along the lines of [3] we observed that Z{k'^) 
became singular at lower but finite k'^. Due to this fact an iterative numerical 
solution proved impossible without constraining the infrared behavior of Z{k^) 
to that according to the discussion of the solution given above. Therefore, we 
conclude that while the expressions for Z{k'^) in [3] with the parameters given 
there are a numerical approximation to the solution of Eq. (16), they do not 
represent the exact analytic result particularly for k"^ <^ ^%cd- 

For the numerical solution of Eq. (26) it turns out to be necessary to study 
the subleading infrared behaviour of F{x) in more detail. This is possible by 
an asymptotic expansion of F{x) at small x analogous to what was introduced 
by Atkinson et al. for Mandelstam's original equation in ref. [2], 

M N 

F{x) = E E a„n^^°+"^°+"+''" . (34) 

m=0 n=0 

The coefficients can be determined from Eq. (26): 

Oon = for n > (35) 

Omo = for m > 

-j^ / m— 1 n— 1 \ 

^mn — I ^m— l,n ^ ] ^ y (^m'n'(^m—m',n—n'C'm—m',n—n' 1 ) ^) 71 > 

with the matrix Cmn given by 

+ 
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2(7o + n7o + n + 2m + 3) 2(70 + n7o + n + 2m + 2) 
9 7 



8(70 + n-'jo + n + 2m + 1) 70 + wjo + n + 2m 

+ ^ 7 ■ (36) 

8(70 + 7170 + n + 2m - 1) 

The double series (34) is thus fixed up to the coefficient aoo of the infrared 
leading term. With this at hand the infrared region up to a matching point 
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xq is treated analytically and the equation for F{x) becomes: 



M N 

EE 



1 + xF{x) n=0 

17 

+ 



7 



2(70 + ^70 + n + 2m + 3) \ x J 
2 9 



Xq 



+ 



2(70 + ^7o + ^ + 2m + 2) V X y 8(70 + 77.70 + n + 2m +1) \ x J 
7 7 



7o + ^7o + 71 + 2m 8(70 + ?T-7o + n + 2m — 1) \ x 



Xq 



a: 



dy f7y^ 17 y'^ 9y 




Similarly the infrared part of the constraint (24) is rewritten with the help of 
the asymptotic series. In the numerical evaluation of the (ultraviolet finite) 
integral in Eq. (24) large momentum contributions are also truncated. Though 
the resulting error becomes negligible for a sufficiently large ultraviolet cut- 
off Xi, convergence is improved by correcting this error analytically. This is 
possible because F{x) can be replaced by its leading perturbative from for 
X > xi and sufficiently large xi, which allows to calculate the finite correction 
term for the integration from Xi to infinity on the r.h.s. of Eq. (24) analyti- 
cally. To obtain this ultraviolet behaviour, the l.h.s. of eqn. (26) is expanded 
for xF{x) 00. The leading terms cancel each other due to the constraint 
(24). The next to leading order terms upon differentiation yield the differential 
equation 



xF'{x) = -7F\x) , 



(38) 



which is solved by 



F{x)^ 



for x — > 00 . 



X 



(39) 



Going to higher orders one starts by expanding (26), neglecting terms of order 
1/x and higher and making use of Eq. (24). Substituting t — Vlnx one arrives 
at 



d 1 r- dgit) 



10 1 

H F^TT7 + 



dtF{t) dt VUt + g{t) 3VUt^ 

where we defined the function g{t) for convenience. Thus 



dgit) _ gjt) ^ 10 1 ^ ^fg'jty 



(40) 



dt 



t 3y/U 



t2 



(41) 
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which is solved by 



c\/l4 10 Int ,^/l\ / 1 > 



3^14 t VtV Ve^^ 

Resubstituting for x finally gives 



V^o Inx («;olnx)3/2 \^(lnx)V2y 



(43) 



where kq = 14 and ni = 70/3. The integration constant c is difficult to deter- 
mine numerically because of the unknown contribution of the terms neglected 
and the extremely large scales necessary for the leading logarithms to domi- 
nate. Wc find roughly that c ^ 1, its exact value is irrelevant, however. With 
X = k'^/X^ and Z{x) = giZ(A;^)/(47r) the leading term in (43) is equivalent to 

Z(e) ^ for e^oo. (44) 

This resembles the renormalization group improved perturbative result, 

Z,t{k') = {Pog'l^kyA'y'^' , (45) 
and allows for the following identifications: 

• The Callan-Symanzik coefficient /3o is given by (3o — kq/ (47r)^, and the above 
value for has to be compared to its perturbative result Kq = ll — {2/3)Nf. 

• Since wc use a momentum subtraction scheme, we identify the QCD scale 
to be used as Amom = A. 

• The anomalous dimension of the gluon field 7a (s') in perturbation theory is 
given by 7a(^) = with 7^ = (13/2-(2/3)A^/)/(47r)^ and the exponent 
in the r.h.s. of (45) is 5/2 — 7°//5o- The fact that 5 = 1 for the asymptotic 
behavior in the Mandelstam approximation confirms that ZgZ^ = 1, which 
is equivalent to 7^ = (3o/2. We thus obtain 7^ = 7/(47r)^. 

The values for the scahng coefficients Po and 7^ obtained from the Man- 
delstam approximation are reasonably close to their perturbative values for 
Nf = 0. The respective values for Mandelstam's original equation (10) are 
Po = (25/2)/(47r)2 and 7^ = (25/4)/(47r)2. 

With the asymptotic infrared as well as ultraviolet contributions integrated 
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analytically, the constraint (24) can now be written, 



M N 
m=0 n=0 



70+n7o+n+2m— 1 




7o + ?^7o + n + 2m - 1 



+ 



dy — ^ + erk \/lnxi 



xo 



y 



14 



(46) 



The integrals in eqn. (37) and (46) are calculated with a Simpson integration 
routine of fourth order. The meshpoints have been chosen equidistant on a 
logarithmic scale, i.e. we have substituted 

j dy — > j duy with u — lny . (47) 

Note that choosing an equidistant mesh in y will not result in a convergent 
iterative process. 

Equation (37) can now be solved iteratively. Starting with a trial function for 
F{x) the r.h.s. of (37) is evaluated and, in order to account for the constraint, 
the result is rescaled appropriately before it is used in the next iteration. To 
leading order the contribution from the infrared expansion is proportional to 
Ooo- It is therefore sufficient to rescale Ooo linearly. This way, the coefficient 
Ooo is determined numerically from the iteration process. For Eq. (12) we 
obtain aoo = 0.29446..., and for (10) its value is aoo = 0.29421.... The infrared 
matching point xq is chosen in a region around 0.2. For this value it proves 
sufficient to truncate the series (34) at M = = 4 (see Fig. 2). For lower 
orders of the series we observe small bumps at the matching point. Changing 
the matching point in a regime from 0.15 to 0.25 has no influence on the 
numerical solution. The lower limit is dictated by convergence of the iterative 
process; the upper limit by a reasonable accuracy of the asymptotic series. 

Using the numerical solution for Z, see Eq. (12), the renormalization group 
invariant product gZ as a function of A;^/A^Qjy[ is shown in Fig. 3. The solution 
to Mandelstam's original equation (10) is almost indistinguishable on this 
scale. 



2.2 The running strong coupling 



With the solution Z{x) the running couphng asifi"^) = Q;s(s) with s = 
/x^/A|^OM follows for the subtraction scheme (15) from Eq. (27). The result 
together with its analytically extracted asymptotic forms is shown in the left 
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0.1 
0.0 

0.0 0.2 0.4 0.6 0.8 1.0 

Fig. 2. The numerical solution F{x) together with several asymptotic series of 
different order. 

0.4 



0.2 

10° 10' 10" 10^ 10^ 

Fig. 3. The numerical solution F{x) and for comparison Z(x) = 1/x + F{x). 

graph of Fig. 4. In the infrared we have, 

as{s) ^47r\ or ^(^.^^ ^ 4^ for /.^ ^ . (48) 

S jJL 

In the ultraviolet one obtains from Eq. (43) the leading logarithmic behavior, 
/ 2x ~ / X 47r / 2Kilnlns 2c \ ( 1 \ ,^fl\ 

This form resembles the two-loop perturbative result if the integration con- 
stant c is set to zero. The value of ki = 70/3 compares to Ki = 51 — (19/3)A^j in 
perturbation theory [16]. The Callan-Symanzik /3 function, (3{g) = ii{dg/dii), 
follows readily for all positive values of the coupling. The numerical result is 





17 



shown in the right graph of Fig. 4. Its hmits are 



Pig) -Po9^ for g ^0 and (3{g) -2g for g ^ oo . 

(50) 

An alternative way to obtain the f3 function follows directly from the cou- 
pling renormalization g^ = Zgg. Because of ZgZ^ = 1 the gluon field renor- 
malization (15) is equivalent to a renormalization condition that relates the 
/^-independent bare coupling to the renormalized one, 

47r Zs An ^ 

— = 47r— = Izifi . 51 

9q 9 9 

This determines the (3 function of the Mandelstam approximation with a mo- 
mentum subtraction scheme (fixing Z{ii^) to a // independent constant as in 
(15)). Since ZgZ^ — 1 the anomalous dimension of the gluon 7^(5') is, 

/X I fi dZs f3{g) 

From the infrared behavior of the renormalization group invariant product 
gDfj_ij{k) we may deduce a linear rising potential between static color sources 
for large distances. This allows us to relate the coefficient of the term in 
the interaction to the string tension a, 

for .^^0. (53, 

With the identification A^qm — (from the ultraviolet) and the infrared 
behavior of the solution, 

gZ{e) ^ 4n^ , (54) 

the non-perturbative result for the gluon propagator relates the string tension 
to the scale Amom, yielding 

2(7 = A^oM (55) 

in the Mandelstam approximation. The string tension can be fixed from quarko- 
nia potentials and Regge phenomenology [17,18]. The result is a value of about 
a = 0.18GeV^. This has also been confirmed in lattice calculations [19]. Here, 
it corresponds to a scale Amom of 600 MeV. 
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10"' 10° 10' 10' 10° 10' 



1 2 3 4 5 

g 



Fig. 4. The strong coupling asd-i'^) with fj? in units of the QCD scale A^jqj^ 
and its asymptotic forms according to (48) and (49) with c = 1 (left); and the 
non-perturbative P function of the Mandelstam approximation (right). 

3 Conclusions 



In this paper we have derived the infrared behavior of the strong coupling 
constant from Mandelstam's approximation [1] to the ghion Dyson-Schwinger 
equation. This approximation results in a gluon propagator diverging like a/k"^ 
for A;^ — > 0. This highly infrared singular behavior of the gluon propagator 
generates a linearly rising potential between static color sources. The corre- 
sponding string tension a can be used to fix the scale in the resulting gluon 
self-energy. In particular, the non-perturbative solution allows to relate the 
string tension to the QCD scale Amom- This relation is in reasonable agree- 
ment with the phenomenological values of the respective physical constants. 
Note that dimensional transmutation occurs here: The coefficient = 2a 
of the infrared singular term is a dimensionful physical quantity, whereas all 
input is dimensionless except for the renormalization scale /x. 

We have shown that the product of the coupling and the gluon propagator, 
gDn„(k), does not acquire multiplicative renormalization in Mandelstam ap- 
proximation {ZgZs = 1). We calculated the gluon self-energy in a renormaliza- 
tion group invariant fashion. Besides the running coupling as{n'^) we obtained 
the non-perturbative Callan-Symanzik (3 function and the anomalous dimen- 
sion of the gluon field for the Mandelstam approximation. In the ultraviolet 
the leading logarithmic decrease of the gluon renormalization function allows 
to identify the scaling coefficients /Sq and 7° at the fix point. The scaling coeffi- 
cients we obtain are in reasonable agreement with the perturbative results. In 
addition, we have shown that there is no further fix point for positive values of 
the coupling in the Mandelstam approximation. The resulting strong coupling 
constant increases without bound in the infrared. 

We have demonstrated how to solve an integral equation with infrared singlu- 
lar solution. This combines numerical and analytic methods based on asymp- 
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totic expansions. It can sucessfully be extended to physically more realistic 
situations in which coupled sets of truncated Dyson-Schwinger equations have 
to be solved allowing for similarly infrared divergent solutions for the propa- 
gators in QCD [4]. 



4 Description of the program 

4.1 The main program 

On startup all relevant variables are initialized. These include the maximum 
number of iterations (NITMAX), the number of meshpoints (MESHPTS) 
used to represent F{x) on a grid, and the orders in the asymptotic expan- 
sion {NMAX, MMAX). Furthermore the infrared matching point a^o and the 
ultraviolet cut-off Xi are set to some appropriate values [XMIN, XMAX re- 
spectively). The variable EPS is introduced to assess the accuracy of the 
result. 

After this, the initial gluon function and the corresponding asymptotic se- 
ries are generated. In the following, in each iteration a new temporary gluon 
function is calculated from equation (37). In order to fullfill the constraint 
(46) this temporary gluon function together with the asymptotic series, i.e., 
the constant aoo, are rescaled appropriately before the result of this two-step 
process is used in the next iteration. 

Convergence is monitored by pointwise comparing the relative deviation be- 
tween the current result for the gluon function and the result of the previous 
iteration. Once the maximum deviation is below the desired accuracy EPS, 
the iteration process is halted. 

Finally, from F{x) the running coupling Q;s(/i^) = g'^{fi'^)/A7i is calculated 
along with the /9-function defined by (3{g) = ii{dg / dii) . The datafile for the 
gluon function (gluon . out) is written in three columns containing (from left) 
X — A;^/A^oM) ^{^) — 1/^ + F{x). The running coupling in the 

format "s = Ai^/A|^oM) <^s(s)", and "g, (3{g)" are written to the datafiles 
alpha. out and beta. out respectively. 

4-2 Subroutines and functions 
Subroutine dijf 

Calculates the first derivative using a five-point formula [20] . 
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Function erfc 

Returns the complementary error function erfc(x) [21]. 
Function Simpson 

Returns the integral of a function which is given at equally spaced abscissas. 
For a sufficient number of abscissas this function uses a closed Simpson rule 
of order l/A^^ [21]. 



5 Testing the program 

Extensive tests were performed to establish resonable ranges to be used for the 
infrared matching point Xq als well as the ultraviolet cut-off Xi. Independence 
of the results on these matching points was obtained within these ranges to 
sufficiently high accuracy. To achieve this, tests have been performed to find 
the necessary order in the asymptotic infrared expansion. The inclusion of even 
higher orders in this expansion was verified to have no considerable effect on 
the results. Similar tests were done to verify the independence of the (suffi- 
ciently large) number of momentum meshpoints as well as the independence 
of the initial values chosen for the gluon renormalization function. 

We solved both, Mandelstam's original equation (10) as well as the improved 
equation (12) obtained by Brown and Pennington in the same truncation 
scheme, with the method presented here. Even though both equations contain 
rather difi^ercnt coefficients, identical qualitative features of their solutions 
were obtained. In particular, the results from the original Mandelstam equa- 
tion were verified to reproduce those of ref. [2] where this equation was solved 
by a different procedure matching the result of an infrared expansion to a 
numerically solved differential equation. 
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TEST RUN 



standard output 

Number of meshpoints : 500 

Order of asymptotic series : M = 4 N = 4 

Infrared matching point : 0.20 

Ultraviolet cut-off : l.OE+08 

eps = l.OE-07 



Convergence achieved after 136 iterations! 
a(0,0) = 2.9446751985E-01 

max. deviation between Fin and Fout : 8.95681E-08 



0.4000000000E-09 
0.4163493887E-09 
0.4333670337E-09 

0.9230051410E+08 
0.9607315655E+08 
. lOOOOOOOOOE+09 



gluon.out 

0.3380217288E-12 
0.3556712162E-12 
0.3742422549E-12 

0.6494579660E-01 
0.6486822745E-01 
0.6479095232E-01 



0.2500000000E+10 
0.2401828914E+10 
0.2307512852E+10 

0.6494580744E-01 
0.6486823786E-01 
0.6479096232E-01 



0.3641128406E-02 
0.3789953965E-02 
0.3944862541E-02 

0.9230051410E+08 
0.9607315655E+08 
O.lOOOOOOOOOE+09 



alpha, out 

. 9478445822E+06 
0.8748652799E+06 
0.8075050070E+06 

0.5300442226E-01 
0.5287788367E-01 
0.5275197512E-01 



. 6327320355E+02 
. 6082727354E+02 
. 5847666993E+02 

0.8161330861E+00 
0.8151583181E+00 
0.8141872450E+00 



beta, out 

-0.1244249422E+03 
-0.1197555950E+03 
-0.1149508614E+03 

-0.4875758004E-01 
-0.4857271919E-01 
-0.4838870479E-01 
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